Polymer translocation into laterally unbounded confined environments 
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Using Langevin dynamics simulations in three dimensions (3D), we investigate the dynamics 
of polymer translocation into the regions between two parallel plane walls with separation R un- 
der a driving force F, respectively. Compared with an unconfined environment, the translocation 
dynamics is greatly changed due to the crowding effect of the partially translocated monomers. 
Translocation time r initially decreases rapidly with increasing R and then saturates for larger R, 
and the confined environment leads to a nonuniversal dependence of r on F. 

PACS numbers: 87.15.A-, 87.15.H- 



I. INTRODUCTION 

The transport of biopolymers through a nanoporc has 
attracted broad interest because it is a challenging prob- 
lem in polymer physics and is also related to many bi- 
ological processes, such as DNA and RNA translocation 
across nuclear pores, protein transport through mem- 
brane channels, and virus injection. Due to its poten- 
tially revolutionary technological applications (T), Q , in- 
cluding rapid DNA sequencing, gene therapy and con- 
trolled drug delivery, a number of recent experimental 
and theoretical (9rl35| studies have been devoted to 
this subject. 

From both the basic physics as well as a technology 
design perspective, an important measure is the scal- 
ing of the average translocation time r with the polymer 
length N, t ~ N a , and the value of the corresponding 
scaling exponent a. Standard equilibrium Kramers anal- 
ysis [36| of diffusion through an entropic barrier yields 
t ~ A^ 2 for unbiased translocation and r ~ N for driven 
translocation (assuming friction to be independent of 
N) [n], HI]. However, as noted by Chuang et al. [ItJ 
the quadratic scaling behavior for unbiased translocation 
cannot be correct for a self-avoiding polymer. The rea- 
son is that the translocation time would be shorter than 
the Rouse equilibration time of a self-avoiding polymer, 
tr ~ N 1+2v , where the Flory exponent v = 0.588 in 
3D and v^d — 0.75 in 2D |37H39| . This observation ren- 
ders the concept of equilibrium entropy and the ensu- 
ing entropic barrier inappropriate for polymer transloca- 
tion dynamics. Chuang et al. [I?) performed numerical 
simulations with Rouse dynamics for a 2D lattice model 
to study the translocation for both phantom and self- 
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avoiding polymers. They decoupled the translocation dy- 
namics from the diffusion dynamics outside the pore by 
imposing the restriction that the first monomer, which 
is initially placed in the pore, is never allowed to escape 
back out of the pore. Their results show that for large N, 
t ~ N 1+2v , which scales in the same manner as the equi- 
libration time but with a considerably larger prcf actor. 
This result was recently corroborated by extensive nu- 
merical simulations based on the Fluctuating Bond (FB) 
[2l| and Langevin Dynamics (LD) models with the bead- 
spring approach l23l |3"H . For driven translocation, Kan- 
tor and Kardar |18| have demonstrated that the assump- 
tion of equilibrium in polymer dynamics in previous work 
[ToL [l2| breaks down more easily and provided a lower 
bound t ~ N l+l> for the translocation time by compari- 
son to the unimpeded motion of the polymer. Using FB 
[22^ and LD [23|,[25[ models, a crossover from r ~ A^ 2 " for 
relatively short polymers to t ~ N 1+v for longer chains 
was found in 2D. In 3D, we find that for faster translo- 
cation processes r ~ A^ 137 [28|, [2j|, while it crosses over 
to t ~ N 1+v for slower translocation, corresponding to 
weak driving forces and/or high friction (30j . 

However, above physical pictures are based on translo- 
cation into an unconfined trans side. Very little atten- 
tion is paid to the dynamics of translocation into con- 
fined environments. It is known that the crowding due 
to macromolecular aggregates and other inclusions in the 
cellular cytoplasm can be as high as 50% by volume and 
has considerable influence on reaction rates, protein fold- 
ing rates, and equilibria in vivo [iol ]. Similarly, polymer 
translocation into a confined environment is subject to 
a large entropic penalty which should dramatically af- 
fect the translocation dynamics. In addition, studies on 
translocation into confined geometries will shed light on 
the dynamics of DNA packaging [4l[ . 

We here investigate the dynamics of polymer translo- 
cation into the region between two parallel plane walls 
(3D) using Langevin dynamics simulations. The paper 
is organized as follows. In section II, we briefly describe 
our model and the simulation technique. In section III, 
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FIG. 1: Schematic representation of polymer translocation 
into a confined environment under an external driving force 
F in the pore. The simulations are carried out in a planar 
confinement (3D), where two plates are separated by a dis- 
tance R. One plate has a pore of length L — a and diameter 

W = 2(7. 

we present our results. Finally, the conclusions and dis- 
cussion are presented in section IV. 



II. MODEL AND METHODS 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains of Lennard- Jones (LJ) 
particles with the Finite Extension Nonlinear Elastic 
(FENE) potential. Excluded volume interaction between 
monomers is modeled by a short range repulsive LJ po- 
tential: U LJ {r) = 4e[(f ) 12 - (f) e ] + e for r < 2 1 / 6 ct and 
for r > 2 1 / 6 cr. Here, a is the diameter of a monomer, 
and e is the depth of the potential. The connectivity 
between neighboring monomers is modeled as a FENE 
spring with Ufene{i") = — ^fci?Qm(l — t 2 /Rq), where r 
is the distance between consecutive monomers, k is the 
spring constant and i?o is the maximally allowed separa- 
tion between connected monomers. 

We consider a geometry as shown in Fig. [TJ where two 
walls with separation R are formed by stationary parti- 
cles within a distance a from each other. One wall has 
a pore of diameter la. Between all monomer-wall parti- 
cle pairs, there exist the same short range repulsive LJ 
interaction as described above. In the Langevin dynam- 
ics simulation, each monomer is subjected to conserva- 
tive, frictional, and random forces, respectively, with [42[ 
mfj = - V{U L j + Ufene) + F xt - fri + Ff , where m is 
the monomer's mass, £ is the friction coefficient, v,; is the 
monomer's velocity, and is the random force which 
satisfies the fluctuation-dissipation theorem. The exter- 
nal force is expressed as F cx t = Fx, where F is the exter- 
nal force strength exerted exclusively on the monomers 
in the pore, and x is a unit vector in the direction along 
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FIG. 2: The radius of gyration at equilibrium state as a func- 
tion of the chain length confined between two walls with sep- 
aration R — 3.5. In the figure, R g , y = R g , z = iZn is the 
direction parallel to the wall, while R g ,x — R± is the direc- 
tion perpendicular to the wall. 



the pore axis. 

In the present work, we use the LJ parameters e and a 
and the monomer mass m to fix the energy, length and 
mass scales respectively. The time scale is then given 
by thj = ima 2 / e) 1 / 2 . The dimensionless parameters in 
our simulations are Rq = 2, k = 7, £ = 0.7 and fc^T = 
1.2 unless otherwise stated. The driving force F is set 
between 1.0 and 15, which corresponds to the range of 
voltages used in the experiments 0, Q . The Langevin 
equation is integrated in time by a method described by 
Ermak and Buckholz [H in 3D. 

Initially, the first monomer of the chain is placed in the 
entrance of the pore, while the remaining monomers are 
undergoing thermal collisions described by the Langevin 
thermostat to obtain an equilibrium configuration. Typ- 
ically, we average our data over 1000 independent runs. 



III. RESULTS AND DISCUSSION 

According to the blob picture [37j , a chain confined be- 
tween two parallel plates with separation a <C R <C R g 
will form a 2D self-avoiding walk consisting of 7i& blobs of 
size R. Each blob contains g = (R/a) 1 ^ monomers and 
the number of blobs is n& = N/g = Nia/R) 1 ^. Thus, 
the blob picture predicts the longitudinal size of the poly- 
mer to be [13, Hi R\\ - n b U2D R - jV«M> <7 (i)'W<'-i ^ 

7V 3/4 (t(^) - 28 , where the Flory exponent v 2 d = 0.75 in 
2D and v = 0.588 in 3D. Fig. [2]shows the radius of gyra- 
tion as a function of the chain length with wall separation 
R = 4.5. We obtain R g>y = R g ^ z - TVO-77+o.oi^ which is 

in good agreement with R\\ ~ TV 3 / 4 . R g x almost does 
not change with N, but R g scales with TV in the same 
way as Rn. The R dependence of R^ is shown in Fig. 
[31 where R g . y — R g . z ~ ^-o.26±o.oi^ m g 00c [ agreement 
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FIG. 3: The radius of gyration at equilibrium state under 
the confinementi between two walls with separation R for 
N = 128. In the figure, R g , y — R g , z = R\\ is in the direc- 
tion parallel to the wall, while R g ,x = R± is in the direction 
perpendicular to the wall. 



with R 



(4) 



0.28 



shows 



For 10.0 > R > 2.5, R g 
linear behavior with R, but R g ~ #-o.26±o.oi_ 

The longitudinal relaxation time tii is denned as the 
time for a polymer moving a distance of order of its 

R 2 

longitudinal size, R\\ . Thus, t\\ scales as t\\ ~ -J- ~ 

N 1+2u2d R 2(1-v 2 dM _ ^2.50^-0.55^ where Q ig the geg _ 

ment length and D ~ 1/N is the diffusion constant. The 
free energy cost of the confined chain in units of fc^T is 
simply the number of the blobs, J- = N(afR) 1 ^ . 

For polymer translocation into confined environments, 
a driving force is necessary to overcome the entropic re- 
pulsion f(R) exerted by already translocated monomers. 
Due to the highly non-equilibrium property of the 
translocation process, it is difficult to estimate the re- 
sisting force f(R). For polymer translocation into the 
region between two parallel plane walls with separation 
R, we assume that the resisting force f(R) scales as 
f(R) = CR 1 for slow translocation processes, with C 
and 7 being the prefactor and the scaling exponent, re- 
spectively. Then, under an external driving force F 
in the pore the translocation time r can be written as 
r ~ F N f(R) ~ fTT^cWTf) w ^ n a being the scaling ex- 
ponent of t with chain length N. Due to ~ for 
the unconfmed system with R ~ oo for slow transloca- 
tion processes, we have 1 — — ~ CR 1 / F. Based on this 
relationship, we can examine the dependence of r on R. 



A. Translocation times as a function of the chain 
length 

Fig. [4] shows the translocation time r as a function of 
the chain length N. For a strong driving force F = 5, 
we find r ~ N a with a = 1.39 ± 0.01 and 1.37 ± 0.01 for 
R = 3.5 and 1.5, respectively, which are quite close to 
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FIG. 4: Translocation time r as a function of the chain length 
N for different R and F in 3D. 
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FIG. 5: The radius of gyration of the chain before transloca- 
tion (t = 0) and at the moment just after the translocation 
(t = r) for F = 15 and R — 3.5. Here, the x direction is 
perpendicular to the wall, while y and z are along the wall. 
The dashed lines are equilibrium values of R g , 7?y and R x for 
chains confined between two walls, as shown in Fig. [2] 



that for an unconfmed system (i? = oo) [29l - l32j . For the 
case with very strong driving force F = 15 and R = 3.5, 
the scaling exponent also does not change, a = 1.37 ± 
0.02. We have checked the translocation velocity with 
respect to the last monomer and found v ~ N@ with (3 ns 
—0.8, which is also the same as that for an unconfmed 
system (3fj|. l3lj. Under a planar confinement, during the 
translocation process the translocated monomers cannot 
have time to diffuse away from the pore corresponding 
to a pronounced non-equilibrium situation. After the 
translocation the distance the last monomer has moved 
is always the radius of the gyration of a chain R g ~ N v a 
for its unconfmed state. Thus, the translocation time 
can be estimated by t ~ — , which is confirmed by our 
numerical results. 
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Above scaling behavior can be understood by measur- 
ing the radius of gyration of the chain before transloca- 
tion and at the moment just after the translocation, as 
shown in Fig. [5j For the chain at the moment just after 
the translocation, R g ~ N - 55 and R g \\ ~ N - 60 , signif- 
icantly lower than their equilibrium values and different 
from the equilibrium scaling of iV 3 / 4 as shown in Fig. 
[3J These results indicate a pronouncedly non-equilibrium 
compression of the chain immediately after translocation, 
as shown in Fig. [5] 

With decreasing the driving force to F = 1.0 and 0.8, 
the scaling exponent remains the same only for longer 
chains. For shorter chains N < 32, the translocation is 
much faster due to considerably lower confinement. 

Decreasing the driving force further to F = 0.5, we 
find the translocation dynamics crosses over to another 
regime with a = 2.44 ± 0.03 at least for N < 64. For 
longer chains, the computation is very expensive due to 
lower translocation probability and a longer transloca- 
tion time. This result is completely different compared 
with the case without the planar confinement (R = oo), 
where we have found r ~ N 1+lJ for F = 0.5 in our pre- 
vious studies [30j . Without the planar confinement, the 
previous prediction by Chuang et al. [TtJ shows that the 
translocation time scales in the same way as the relax- 
ation time of the chain for unbiased translocation. Here, 
our exponent a = 2.44 ± 0.03 for F = 0.5 is very close 
to 2.50, which demonstrates that r scales in the same 
way as the 2D longitudinal relaxation time for the weak 
driving force. The reason is that the weak driving force 
F is almost balanced by the resisting force due to the 
confinement entropy and the net driving force is close to 
zero. Compared with strong driving cases, the chain con- 
formation is not so compressed, as shown in Fig. [6] Here, 
we also should point out that the exponent a = 2.44 is 
a little smaller than 2.50 due to the fact that during the 
translocation process the untranslocated monomers are 
not confined. To obtain the exponent 2.5, one should 
initially squeeze the whole chain into the confined spaces 
between two plates and then record the time during the 
chain moving a distance of order of Ru . 




10 15 20 



FIG. 6: The chain conformation projected onto the yz plane 
for N = 128 under strong and weak driving forces, respec- 
tively. 



that before the translocation. At the moment after the 
translocation, R gtV — R g , z ~ R~ 20 for R < 9, which is 
slightly lower than the equilibrium exponent —0.28. In 
addition R g . x ~ i? 1 ' 43 for R < 3, somewhat lower than 
the equilibrium exponent 1/v = 1.70. For larger R, such 



as R = 19, R g 



Ra 



R a 



B. Translocation times as a function of R 

Fig. [7] shows the R dependence of the translocation 
time for F = 5.0 and N = 128. Initially, r decreases 
rapidly with increasing R and then almost saturates for 
larger R. Quantitatively, we check 1 — — MS ir.l function 
of R, see the insert of Fig. [7J For R < 3, we find 7 = 
— 1.61 ± 0.06. It is close to the exponent — = —1.70, 
which is the theoretical value for equilibrium transloca- 
tion processes with f(R) ~ R~ x l v . This result also shows 
that for R < 3, the chains really feel strong confinement 
and the translocation is slow enough for the chain to be 
close to local equilibrium. Fig. [5] shows the radius of gy- 
ration of the chain as a function of R, where R g ~ R~ 20 
at the moment after the translocation is still lower than 



C. Translocation times as a function of the driving 

force 

According to the Langevin equation, when the iner- 
tia term can be neglected, the balance of the frictional 
force and the driving force, £u = F, gives r ~ £/F for 
the case without planar confinements. The Langevin dy- 
namics simulation results [23l l30j confirms r ~ 1 j F only 
for slow translocation with lower values of F/£. For fast 
translocation with higher values of F/£, we have found 
t ~ F~ ' 8 due to the highly deformed chain confor- 
mation on the trans side, reflecting a pronounced non- 
equilibrium situation [3(|. Fig. |H] shows the plot of r 
as function of F for translocation into a region of a pla- 
nar confinement. Under strong driving forces, such as 
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FIG. 7: Translocation time r as a function of R for chain 
length N = 128 under F = 5 in 3D. 
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FIG. 8: The radius of gyration of the chain before transloca- 
tion (t = 0) and at the moment just after the translocation 
(t = r) for N = 128, F — 5 and different R. Here, the x 
direction is perpendicular to the wall, and y and z are along 
the wall. The dashed lines are equilibrium values of R g , R\\ 



and R x 
Fig.E] 



for chains confined between two walls, as shown in 



F > 5, compared with the case without confinement 
r ~ ^-o.80±o.02 j g s ^jjj observed. Under a planar con- 
finement, it is more difficult for translocated monomers 
to diffuse away from the pore, which slows down the 
translocation. As a result, the exponent must be lower 

than — 1. Under weak driving forces F < 2.5, we find 
r <~ j?-i-25±o.03_ j a ^ a can ^ e g ttc j w j tn an em pi r _ 

ical function r = 2300F- 125 [1 + (f/3) 10 ] 045 , yielding 
the details of the crossover from the weak force to the 
fast strong force regime for this case. The seemingly high 
exponent 10 in the empirical function is necessary to ac- 
count for the relatively fast turnover between the scaling 
t ~ F -1 - 25 at weak forces and the behavior t ~ F~ ' 8 at 
strong forces, as visible in the shown data. For translo- 
cation into two parallel plates with separation i?, the 
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FIG. 9: Translocation time r as a function of F for N = 128 
and R — 3.5 in 3D. 



effective driving force is F(l — Cf(R)/F). For weaker 
driving forces, the resisting force f(R) becomes more and 
more important because (1 — Cf(R)/F) is closer to zero. 
Therefore, with decreasing F the translocation is greatly 
slowed down, resulting in stronger dependence of r on F. 
In our previous results (28j , we find that the polymer-pore 
interactions can also lead to this behavior. 



D. Waiting time distribution 



The time of an individual segment passing through the 
pore during translocation is an important quantity con- 
siderably affected by the translocation dynamics. The 
noncquilibrium nature of translocation has a significant 
effect on it. We have numerically calculated the waiting 
times for all monomers in a chain of length N. We define 
the waiting time of monomer s as the average time be- 
tween the events that monomer s and monomer s + 1 exit 
the pore. Fig fTUf a) shows the waiting time distributions 
for different chain lengths for R = 3.5 and F = 5. One 
obvious feature is that before reaching their maxima the 
waiting times almost follow the same pathway, which is 
dominated by the strong driving force. In addition, the 
maxima occur at s max > A/2, such as s max ~ 0.65-/V 
and 0.677V for N = 96 and 128, respectively. Fig. ITOTb) 
shows the waiting time distributions for N = 128 and 
F = 5 for different R. As expected, with decreasing 
R, the waiting times increase. Fig llll shows the wait- 
ing time distributions for A = 64 and R = 3.5 under the 
weak driving force F = 0.5. After half of the polymer has 
been translocated, it still takes a considerably long time 
to exit the pore for the remaining monomers, compared 
with strong driving cases shown in Fig llOl This implies 
different translocation dynamics in the two regimes. 
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IV. CONCLUSIONS 



Using Langcvin dynamics simulations, we investigate 
the dynamics of polymer translocation into the region be- 
tween two parallel plane walls with separation R. Com- 
pared with the chain passage into an unconfined envi- 
ronment, the translocation dynamics is greatly changed. 
In particular, the translocation time r initially decreases 
rapidly with R and then saturates for larger R in 3D, and 
the confined environment leads to nonuniversal depen- 
dence of r as a function of the driving force. For polymer 
translocation into 3D confinements with R = 3.5, we find 
that under the weak driving force F = 0.5 the translo- 
cation time r scales with chain length N as r ~ N 2Ai , 
which is in the same manner as the relaxation time of a 
chain under a planar confinement, while it crosses over 
to N as r ~ TV 1 ' 39 for fast translocation. These be- 
haviors are interpreted by the waiting time of individual 
monomer. 
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FIG. 10: Waiting time distribution for different (a) N and (b) 
R. 
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FIG. 11: Waiting time distribution for F = 0.5 and N = 64. 
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